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We consider integer quantum Hall states and calculate bulk entanglement spectrum by formulating 
the correlation matrix in guiding center representation. Our analytical approach is based on the 
projection operator with redefining the inner product of states in Hilbert space to take care of 
the restriction imposed by the (rectangle-tiled) checkerboard partition. The resultant correlation 
matrix contains the coupling constants between states of different guiding centers parameterized by 
magnetic length and the period of partition. We hnd various band-crossings by tuning the flux <f> 
threading each chekerborad pixel and by changing filling factor u. When v = 1 and <f> = 2tt, or 
V = 2 and $ = tt, one Dirac band crossing is found. For u = 1 and $ = tt, the band crossings are 
in the form of nodal line, enclosing the Brillouin zone. As for v = 2 and $ = 27r, the doubled Dirac 
point, or the quadratic point, is seen. Besides, we infer that the quadratic point is protected by C 4 
symmetry of the checkerboard partition since it evolves into two separate Dirac points when the 
symmetry is lowered to C 2 . In addition, we also identify the emerging symmetries responsible for 
the symmetric bulk entanglement spectra, which are absent in the underlying quantum Hall states. 

PACS numbers: 73.43.Cd,03.67.Mn,73.43.Nq 


I. INTRODUCTION 

Entanglement spectrum, which provides finer charac¬ 
teristics for the many-body entanglement than the en¬ 
tanglement entropy, has been used to characterize the 
topological quantum states.^^ One characteristic for these 
states is the existence of the gapless mode in the entan¬ 
glement spectrum.l^H^ This reflects the entangled nature 
of the underlying quantum state as there is no dominant 
term in the wavefunction to be approximated as a prod¬ 
uct state. On the other hand, gapless excitations at edges 
are the main character distinguishing the gapped systems 
that are topological from those that are not. The Dirac 
band crossings found on the surface of topological insu- 
latorJ^ are examples of electron systems. The partition 
boundary for evaluating the entanglement spectrum, al¬ 
beit artificial, plays the similar role as the physical edge 
hosting the topological gapless edge modes. Then, the 
correlation between entanglement and the topological or¬ 
der motivated the worlP to demonstrate a general rela¬ 
tion between entanglement spectrum and the physical 
edge state spectrum of topological states. 

The entanglement spectrum, by construction, depends 
on how we partition the ground state in either real space 
or internal space. The simplest partition is just to cut 
the system into two halves by a single boundary in a 
real space, which is also adopted in most of the studies 
mentioned above. This partition usually breaks some of 
the underlying translational or discrete symmetries, es¬ 
pecially for the systems living on the lattice. It is then 
interesting to see what the new feature of the entangle¬ 
ment Hamiltonian and the relation with the physical edge 
Hamiltonian of the topological state will arise if we par¬ 
tition in a different way to preserve some lattice struc¬ 
ture. Because of the residual lattice structure left by the 
partitioning, one naturally expects the emergence of the 


band-like structures for the entanglement spectrum, and 
the ensured gapless modes for the entangled topological 
state lead to the expectation of the band-crossings. 

Along the line, the worlJi^ studied the integer quan¬ 
tum Hall ground state on a square lattice and its entan¬ 
glement spectrum under a checkerboard partition. By 
doing so the translational symmetry is reduced to that 
of a superlattice, and their remarkable findings include 
a family of entanglement spectra which shows how the 
gaps are evolved to zero at a point when the equal and 
extensive partition is employed. Under such extensive 
partition, the bulk entanglement spectrum (BES) is thus 
coined. The subsequent worlP^ has also confirmed the 
the existence of a Dirac band crossing at T point when 
total Chern number of the underlying ground state is one. 
The debate on whether a gap is open or a quadratic band 
crossing emerges when the total Chern number is even 
can be clarified if further analytical investigations are 
taken, which is one purpose of the present paper. More 
recently, implementation of extensive partition to one¬ 
dimensional interacting system has also been made.l^^lti^ 

The topological insulators/superconductors are typi¬ 
cally characterized by the patterns of the Dirac band 
crossings at the physical edges, the robustness of which is 
protected by some discrete symmetries s uch as particle- 
hole, time-reversal and chiral symmetriesBased on 
the similarity between the entanglement spectrum and 
physical edge spectrum for the topological states, one 
may expect that the band-crossings of BES could be also 
protected by some discrete symmetries. It is then in¬ 
teresting to find out the explicit form of entanglement 
spectrum for some underlying topological quantum state 
and then identify the emerging discrete symmetries. 

In this paper, we address the above issues by consider¬ 
ing the integer quantum Hall states under a checkerboard 
partition in continuous space, which allows an analyti- 
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cal investigation and gains more insights than previous 
numerical works do. Given the exact relation between 
correlation matrix for free-electron system and entangle¬ 
ment spectrumfi^H^ it suffices to focus on the correlation 
matrix. Hence, we introduce the notion of guiding cen¬ 
ter representation in which the inner product of states is 
redefined to take care of the underlying partition. The re¬ 
sultant correlation matrix resembles the one-dimensional 
hopping Hamiltonian. Besides, the matrix is always in 
the form of an identity matrix times ^ plus a purely imag¬ 
inary matrix, from which we immediately conclude that 
the complex conjugation is the particle-hole symmetry of 
BES. We further demonstrate the BES for v = \ and 
V = 2 cases and find the Dirac and quadratic band cross¬ 
ings when <i> = 27r, which is similar to the spectra in single 
layer graphen^^ and bilayer graphene.!^ In addition, the 
quadratic point is inferred to be protected by C 4 symme¬ 
try since it is evolved into two separate Dirac points as we 
lower the symmetry of partition into C 2 . The breaking of 
one quadratic point into two Dirac points also occurs in 
the nematic ordered states in bilayer graphene.^ More¬ 
over, the lone Dirac point can also be found in ^ = 2 and 
$ = TT, while the band crossing become the nodal line in 
the case of v = 1 and $ = tt, which is reminiscent of the 
nodal structure in Weyl semimetal.^ 

The paper has the following organizations. In Sec. 
the notion of guiding center representation for correla¬ 
tion matrix is introduced with a simpler example in left- 
right partition. The formulation to obtain the correla¬ 
tion matrix in checkerboard partition then follows. In 
Sec. |HI[ the flux through each square is taken as 27r and 
we demonstrate the BES of ^ = 1 and v = 2 integer 
quantum Hall states where the band crossing points are 
found. The associated discrete symmetries of the entan¬ 
glement Hamiltonian are also identified. We also consider 
the case of half-flux and obtain the nodal lines in BES. 
In the end, the conclusions is drawn in Sec. [V] 


II. CORRELATION MATRIX IN GUIDING 
GENTER REPRESENTATION 

To obtain the entanglement spectrum for a underlying 
quantum state I'k), we should first construct the reduced 
density matrix for the subsystem A. It is given by 
PA '■= TrB|4')('I'| where we have traced out the Hilbert 
space of the subsystem B if we partition the total sys¬ 
tem into subsystems A and B. Then the entanglement 
Hamiltonian is := — lop a whose eigenvalues form the 
entanglement spectrum. 

The entanglement Hamiltonian is usually nonlocal and 
hard to obtain in closed form even for gaussian systems. 
However, in the case of free fermion system, the usual 
approach to obtain the entanglement spectrum, the set of 
eigenvalues of entanglement Hamiltonian is through 
the correlation matrix C with the aid of the following 
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FIG. 1. (color online) The checkerboard partition with each 
pixel of side length a is employed in the paper to study the 
“band” structnre of bulk entanglement spectrum. 

equalitjEHSl 

= ( 1 ) 

in which denotes the transpose of He- The informa¬ 
tion of partition of system is clearly seen in the coordinate 
representation for correlation matrix C'(r,r') with r and 
r' both lie within the focused subsystem, say A. In fact 
C'(r,r') is the matrix ele ment of the projection operator 
onto the occupied states,l^ ^^ * ^^ l 

p= l“)H. (2) 

E^<Ef 

in which the eigenstate |a) of physical Hamiltonian 
is indexed by the quantum number a. In a lattice, 
r and r' shall run through a subset of sites, so the 
correlation matrix C is a finite one. In continuous space, 
one may write down C in coordinate representation. 
The eigenvalue problem of C then becomes an integral 
equation. A simple example is available from the integer 
quantum Hall stat^^ with left-right partition. It 
was showrP that for filling factor 1 / = 1 the correla¬ 
tion matrix in coordinate representation is C(r, r') = 
exp [-{y - y'f/^ - {x- x'Yj^ - i{x -f x’){y - y')/2]/27r, 
and the following integral equation can be solved ana¬ 
lytically. 


/ dr'C{r,r')f{r') = af{r) , (3) 

J r,r'G A 

with A denote the subsystem. In fact, the eigenfunction 
/of Eq.[3] coincides with the lowest Landau level wave- 
function, and the corresponding eigenvalue a is equal to 
the probability of locating / in the subsystem A. 

The approach of solving the integral equation such as 
Eq .0 becomes much more difficult when the partition is 
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not the left-right one. Now we shall demonstrate that 
Eq-i is equivalent to the following eigenvalue equation 
of projection operator P, 


P\a) = a\a) , (4) 

along with the particular dehnition of inner product of 
states in Hilbert space, 

((a|a')) := / drf*{r)fa'{r) , (5) 

JreA 

associated with any pair of eigenstates |a) and |a') of 
physical Hamiltonian. The above dehnition of inner 
product differs from the original one in that the range 
is restricted to the subsystem A of partition. Without 
the restriction, the eigenstates of physical Hamiltonian 
are also eigenstates of P with eigenvalue being either 1 
or 0 depending on wether it is occupied or not.^ How¬ 
ever, it is the restriction that makes the complete set 
{|a)} not an orthonormal one, and hence the eigenvalues 
are generally 0 < a < 1. 

Now we introduce the notion of the guiding center rep¬ 
resentation for correlation matrix. Assuming the size of 
the two-dimensional space is x Ly, the Landau level 
states are represented by the ket \n,X) with energy in¬ 
dex n = 0,1,2,... and X being guiding center coordi¬ 
nate. If the Landau gauge is used, the wavefunction is 
^'nx(r) = {l/Lyy^^exp{iXy/^)(j)n{x-X) where 4>n{x) 
is a one-dimensional harmonic-oscillator eigenstate and 
£ = (hcjePlY^"^ is the magnetic length. The projection 
operator associated with hlling factor v reads. 


\nX){nX\ . (6) 

n=0 X 

First consider v = 1 and the left-right partition. It 
is clear that ((0,A'|0,A)) oc 5x,x' as a result of the 


preserved translational symmetry along the boundary 
of partition which is chosen to coincide with the y di¬ 
rection in the Landau gauge. Thus, applying Pi to 
|0,A) gives the eigenvalue in terms of error function 
dx = ((0,^|0,A)) = [1 -I- Erf(A/.^^)]/2, which is pos¬ 
itive and less than unity. In particular, uq = 1/2 shall, 
according to Eq. correspond to a zero in entanglement 
spectrum. Hence the approach of Eq. [^is a valid one. 

For a more general partition, the eigenvalue equation 
Eq. following from Eq. becomes 


E r^ (a) (cn) 

Cn,X-n',X' S^, x' — Sn,X > 

n',X' 


( 7 ) 


in which s denotes eigenvector in the guiding center rep¬ 
resentation, i.e., := ((n,A|a)), and the correlation 

matrix in this representation is 


:= ((n,A|n',A')) . (8) 

Here we have adopted the inner product defined in Eq. 
so that C is not the identity matrix. Consequently, a is 
simply the eigenvalue of C in Eq. and the correspond¬ 
ing entanglement energy is given by = In . 

With the general approach outlined above, we are 
ready to focus on the entanglement spectrum associated 
with the checkerboard partition shown in Eig. Unlike 
the left-right partition in which the matrix C is diago¬ 
nal with respect to guiding center A, the checkerboard 
partition allows the couplings between different guiding 
enters X and X' as long as the separation between them 
is an odd multiple of G£‘^ with G = xja the reciprocal 
lattice vector. The square pixel has the area of a^. This 
can be understood from the explicit calculations of C, 
which reads 


Gn,X-n',X' = j dv[Ui{x)Ui{y) + U2{x)U2{y)]^nx{^)^n'X'{x) 


i. . , 5x',X+(2m-l)GP 

= A^n,n'Ox,x' + * ^ -—TX-A; n , A ) , 

/ ( }m — 1 ITT 




(2m — l)7r 


(9) 


in which the restriction to the region in white (Fig. is replaced by the series of products of Heaviside functions 
Ui(x) = X)mGZ { 0 ( 3 ; — 2ma) — 0[a: — {2m + l)a]} and U 2 {x) = Ui{x + a). The diagonal term of 7 appears due to the 
fact that all white squares collectively occupy half of the total space and that the probability p is uniform along 
y direction. The constant A appeared in the second line is a real number, which can be calculated with the following 
integral. 


A(n, A;n', A') = J dx [C/i(a;) — C/ 2 (x)](/„(a; — A)(/„/(a; — A') 


E 


{n’)\ ^ piT 


exp 


ipG 


X + X' 




pG, 


A'-A 


£2 


( 10 ) 

( 11 ) 
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in which the symbol S refers to taking the imagi¬ 
nary part of what follows. Thus, the correlation ma¬ 
trix C in guiding center representation resembles a one¬ 
dimensional hopping Hamiltonian. We may expand the 
square wave into the Fourier series, namely Ui{x) — 
U 2 {x) = X]p[4sin(pGa::)]/(p7r) over p = 1,3,5,.... The 
overlap of two separate Gaussian functions is often en¬ 
countered in the context of two-dimensional interact¬ 
ing electr ons in the presence of perpendicular magnetic 
f^eidpZUll and the function = {n\/n'\)[(—qy + 

ig,,)£/y2]’"'-"exp(-g2£2/4)/:;j'-"(g2^V2) with C the 
associated Laguerre polynomial. 

It is of crucial importance to note that matrix C in 
guiding center representation is Hermitian and that the 
off diagonal part G — is pure imaginary. Therefore, 
the entanglement Hamiltonian TJ' = ln[(I — G)G“^] has 
the particle-hole symmetry under /C(I — G)/C = G with 
/C the complex conjugation. 


0 G/2 2a 4a 

FIG. 2. The continuous space of guiding center X can be 
measured by the two length scales 2a and 

defined within the first Brillouin zone [0, G] x [0, G]. With 
this basis labelled by k, the G matrix is simplified into a 
{vN) X {vN) one, leading to vN bands in the BES subse¬ 
quently. In the followings, we shall specifically calculate 
the cases for v = 1 and v = 2 and study the correspond¬ 
ing BES. 


A. u=l 


III. ENTANGLEMENT SPECTRUM 

Now we are in the position to obtain the entangle¬ 
ment spectrum from the matrix G in Eq. Essentially, 
G is a huge matrix and can be thought of as motion 
of hopping in the one-dimensional space of X shown in 
Fig.i Moreover, as a consequence of the checkerboard 
partition, hopping only occurs between X and X ± G£^, 
X ± 3G£^, and so on. In addition, it is obvious that G is 
invariant under the following action of shift by 2a, 


{{n,X\n', X')) = {{n, X + 2a\n', X' -b 2a)) , (12) 

which can be seen from the sine function in Eq. as 
well. This periodicity simplifies the the problem since we 
can change to the Bloch basis. Now, as seen in Fig. 
we may parameterize the guiding center coordinate as 

X{xo,iJ) = xq+ 2ia +, (13) 

with 0 < a;o < G£^ and i gZ. The second integer j G Z^v 
with 


N := 


2a 



(14) 


which is determined by the flux 4) = a^/£^ threading 
each pixel in Fig.[^ It is clear that the index j labels the 
X bands. By Fourier transforming with respect to the 
2a-periodicity, we obtain the following Bloch basis. 


|n,A) hG |n,k,j) = J—J^e^'^‘^^ln,X(ky£^,l,j)} , 

' ^ lez 

(15) 

which shall acquire a phase factor exp (ikx2a) upon trans¬ 
lating one period along the x axis. The momentum k is 


We consider the situation that the flux threading each 
pixel is 27r, which corresponds to A = 4 according to 
Eq. |14| In this case, we may suppress the energy index 
n and the matrix elements ((k', j'|k, j)) = i5k,k'Gk. The 
orthogonality follows from the Bloch basis in Eq. and 
the relation in Eq. Eor convenience, we may rearrange 
the column vector such that the order of j is (0, 2,1,3). 
It can be shown that the 4x4 matrix is 


/i 0 A B\ 
^ _ 0 i Q i? 

A* Q* I 0 
\b* R* 0 i / 

and the matrix elements read. 


(16) 


A = ^((X( V", 0, 0)lX(ky£^, I, l)))e*'=^(4'+FG«'(I7) 

I 

B = ^((X(V^0,0)|A(fcy£2,;,3)))e*'=='(^'+3)^^'(I8) 

I 

Q = ^((X(V^0,2)|A(V^^,l)))e*"^^"'■'^''''(19) 

I 

R = '^{{X{ky£^, 0, 2)\X{ky£'^, /, 3)))e*'="(^'+F‘=^'(20) 


Note that the associated phase in above sums is deter¬ 
mined by the separations between guiding center coor¬ 
dinates. In obtaining the above relations, we have used 
Eq. 12 and the fact that 2a = 4G£^. In addition, the 
coefficients are not all independent, and it can be shown 
that R = —A and Q = —B based on the observation 
that ((A(fcy£^,0,j)|A(fcy£^,^,/))) changes sign if j and 
{41 + j') are simultaneously shifted by 2. 

After some algebra, the two independent coefficients 
for general k read 
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FIG. 3. (color online) The bands in entanglement spectra 
i5e(k) of integer quantum Hall state are shown for the parts 
close to zero energy. Dirac point and quadratic point appear 
at the origin F point in the case of iz = 1 (left) and v = 2 
(right), respectively. 


are the degenerate eigenstates of Ck=o with eigenvalue 
1/2. The two-component spinors \x, ±) are eigenstates of 
with ± eigenvalue, respectively. One may immediately 
see that the two states ipi ^2 correspond to the zero modes 
of He through Eq. 


B. = 2 


Following the above approach, the double-bared Ck 
denotes the corresponding 8 x 8 matrix for v = 2 integer 


quantum Hall state in the basis of Eq. 15 It has the 
form. 






pqiT^ 


( 21 ) 


^ = - 4 *E 


Q 177 Q kx 




P,(l 


pqiT^ 


,( 22 ) 


in which both sums are over p = 1,3,5,... and q = 
1, —3, 5, —7,.... Note that the rescaled momentum k := 
k/G is dehned. Besides, we may also write 


Gk = ^h+<^{A)Ni+Q{B)N2+^{A)Mi+^iB)M2 (23) 

with the representation Ni = —ay 0 Uz, N 2 = —a^ 0 ay, 
Ml = ax®az, and M 2 = —ay®ay. In order to write down 
the low-lying part of BES analytically, we only retain the 
largest terms, which correspond to p = q = 1, in A and 
B. Eigenvalues of Gk are shown to be. 


Gk = 


Goo Goi 
Gio Gn 


(27) 


which consists of the 4x4 matrices, Goo and Gn, the intra 
Landau level couplings, and Goi and Gio ,the inter Lan¬ 
dau level couplings. Gqo is the same as Gk in Eq. IT The 
difference between Gn and Goo is that additional factor 
due to Laguerre polynomial G? enters the off-diagonal 
part of Gn. It is clear that pi ^2 are also eigenstates of 
Gn. Therefore, in the absence of inter-Landau level cou¬ 
plings Goi, two Dirac cones shall appear in BES at T 
point with the ratio of respective “Fermi velocity” to be 
rj := Ci{G‘^i‘^) = {tt/2 — 1 ), approximately. 

Now we are going to show that the BES associated 
with Gk has a quadratic band crossing point at k = 0 . 
The situation here is very similar to the bilayer graphene 
in the sense that the way the two Dirac cones are merged 
is determined by the “interlayer” coupling Goi, which has 
the general form. 


Ok 


1 4 

= i±e,(k) 


sin^ 0+ -I- sin^ 0_ ± 2 sin 9+ sin cos( 7 rfca;) 


1/2 


(24) 


/ 0 Q A B' \ 

0 0 Q' B' ] 

A" B" 0 0 

\Q" R” 0 0 / 


(28) 


which gives rise to four bands that are symmetric with 
respect to The angles are 9± = ± |). Then BES 

is obtained from Ee(k) = ln[(l —Q;k)/ak]- Fig. [^displays 
the numerical results of i5e(k). The Dirac band crossing 
at zero energy results from taking the plus sign in the 
bracket, namely Eg « ±4e+(k) « ±8-\/2 exp(— 7 r/ 4 )fc/ 7 r. 

The momentum k = 0 is a special point at which A = 
B and = ^B = 0. We may thus write 


Gk=o — 2^4 + 


fjy(l2 fXx) 


-ayih + ax) 


with which we may show that the pair of states, 


(25) 


Unlike Goo, Gqi has vanishing diagonal elements be¬ 
cause of Eq. 1^ Besides, the essence of the presence 
of a quadratic band crossing point is that at k = 0 
all the primed coefficients A, B', Q', and R' van¬ 
ish while the double-primed ones do not. It could be 
understood by hrst inspecting A, which is a sum of 
((n = 0,X(0, 0,0)|n = 1,X(0, Z,l))) over integer 1. On 
the other hand, the coefficient A' is a sum of ((n = 
0,X(0, 0, l)|n = 1, X(0,/, 0))). Based on Eq. 11 it can 
be shown that at k = 0 


A {A') cx E 

P,Q 


pcospqj T qsinpqj 
pq ® 


(29) 


ipi = 




in which the sum is over p = 1,3,5,... and q = 
(26) 1, —3, 5, —7,... as before. Note that the minus (plus) sign 

is used exclusively for A {A"). The appearance of both 
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cosine and sine in above expression is due to the fact that 
the F function in Eq. for inter-Landau level coupling 
is complex. Now one can easily see that when p = \q\ the 
contributions vanish in A. In addition, the cross terms, 
for example the pair of contributions from (p, q) = (1, —3) 
and {p,q) = (3,1), are exactly cancelled in A'. In con¬ 
trast, the coefficient A' is not zero at k = 0. It can be 
shown that B' vanishes at k = 0 in the same manner. 
Now we only retain the dominant terms in each coeffi¬ 
cient, 

A' = -4isin(^ky)exp(i^k^)-^^ = -(B')* , (30) 

and 



FIG. 4. (color online) Upper panel: BES of = 1 and <1? = tt. 
Two lines of node appear at fea, = ±y and ky = ±^, and 
the BES respects the C 4 symmetry of checkerboard partition. 
Lower panel: BES of = 2 and <1> = tt. This is equivalent 
to = 1 and = 27r, and the appearance of a single Dirac 
band crossing located at (^, y) confirms the manifestation of 
Chern number in the band crossing in BES. 


A" = _4*cos(^fc,)exp(-*^fc,,)^ = -{By . (31) 


D. Spectrum symmetries 


Plugging A'(k = 0) = S"(k = 0) « t = -7=16-7 into 

Ck, the BES oi 1 / = 2 state is obtained and is shown in 
Fig. [3} Moreover, one can show with some algebra that 
the low-lying part of BES is approximately the solution 
of. 


det 


2t 

4 Ee 
2t* 

4 'I 


= 0 . 


(32) 


and it is reminded that e_|_ oc fc gives rise to the Dirac 
point in 14 = 1 case. Here when |t| ^ e+ the resultant 

quadratic band crossing in BES is ^2^ 


C. TT-fluX 


Here we consider $ = tt, which corresponds to a ?= 
TT and from Eq. 14 we obtain N = 2. Thus the corre¬ 
sponding correlation matrix C is a 2x2 one. For J4 = 1, 
it can be shown that C = ^^2 + Aa+ + H*cr_ with 


A = 4i 


sin[p7r(fcy 






[p^W] 


PQTT^ 


(33) 


in which the sum is over p = 1, 3, 5,... and q — ±1, ±3,.... 
It can be shown that H = 0 along the lines kx = ^ and 
ky = ^. Therefore, as shown in Fig.|^ two lines of node 
appear in the corresponding BES which respects the C4 
symmetry. 

Similar calculations can be carried out for the case of 
V = 2 and $ = tt, which is topologically equivalent to 
the case of i4 = 1 and $ = 27r. The equivalence is further 
confirmed by the presence of a single Dirac crossing at 
k = ^(1,1) shown in Fig. 


As shown in Fig. both BES’s are symmetric with 
respect to zero. Here we shall explore the corresponding 
discrete symmetries in analogy with the time-reversal, 
chiral, and particle-h ole sy mmetries in the mean-field 
fermionic Hamiltonian.E^E^ We start with v = 4 case. 
For $ = 27r, let us focus on the matrix Ck in Eq. 
First examine the action of A. Since lC{i)K, = —i and 
A(k)A = — k, one can see that A in Eq. and B in 
Eq. are invariant under complex conjugation with k 
defined within BZ. Therefore, the antiunitary operator 
It = 1C commute with Ck and we may regard it as the cor¬ 
responding “time-reversal symmetry” in this representa¬ 
tion. Moreover, the matrix Ck — has only off-diagonal 
blocks. Therefore, E = Uj 0 I2 is the “chiral symme¬ 
try” of BES since E-^CkS = 1 — C^. Lastly, it is clear 
that the product S = E/j shall serve as the particle-hole 
symmetry operator. 


23 


For 1 / = 2, the the diagonal blocks, Coo and Cn, of the 
matrix in Eq. |27| share the same spectrum symmetries of 
Ck in Eq. |23[ So the above three symmetry operators 
are also the spectrum symmetries of Cqo and C\\. Fortu¬ 
nately, the off-diagonal block Coi has vanishing diagonal 
elements, so it also anticommute with E. Consequently, 
the chiral symmetry is accounted for with the operator 
I2 ® E. However, since It acts differently on the coeffi¬ 
cients A in Eq. 30 and A' in Eq. 31 the time-reversal 
symmetry is broken. 

As for 14 = 1 and 4) = tt, it can be seen that the matrix 
C — |l2 commutes with /* and anticommutes with the 
chiral symmetry operator The product a^It reflects 
the BES. 


Whether the band crossings in BES are protected can 
be demonstrated by adding perturbations that break 
the above symmetries to, for example, Eq. fol¬ 
lowing the framework of c lassifying topological insu- 
lators/superconductors.*^^*^ However, one has to keep 
in mind that these additions may not correspond to 
tunnable parameters in physical system, and only for the 
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g=l.05 



FIG. 5. (color online) BES oi v = 2 and $ = 27r with con¬ 
sideration of anisotropic factor g ~ ^ = 1.05 where a and 
b are length of sides in x and y directions, respectively. The 
quadratic band crossing point appeared when g = 1 is now 
split into two Dirac points. 


purpose of testing the robustness of band crossings. 


IV. ANISOTROPIC PARTITION 


When the pixel in Fig. is of rectangle shape with a 
ration of two sides g = ^ > 1, the two reciprocal vec¬ 
tors Gx = 'Kja and Gy = 7r/6 are not equal. It follows 
that {GxG)^ = gT^^ 1^ and (Gyi)^ = Thus, the 

product GxGyi^ = 7r/2 for the flux d) = 27r remains. The 
breaking of C4 symmetry into a lower C2 symmetry does 
not ruin the symmetry of intra-Landau level coupling 
constants. For ^ = 1, it can be seen from the coefficients 


in Eq. 21 and 22 that the condition A — B still holds at 


k = 0 despite the fact that the factor of exp [— | -|- q'^)] 

in the sum is changed into exp [—f (pp^ + g~^q^)]- Con¬ 
sequently, the Dirac band crossing in BES is intact, inde¬ 
pendent of the anisotropic factor g. The low-lying spec¬ 
trum is still isotropic in terms of the rescale momenta 
kx = kx/Gx and ky = ky/Gy. 

It is interesting to consider the case oi v = 2 and ex¬ 
plore the effect of g on the quadratic band crossing. Pre¬ 
viously, the vanishing of coupling constants A! and B' at 
k = 0 is essential for the emergence of quadratic point. 
However, it is no longer the case when g ^ 1 since the 
coupling constants A' and B' are, to the lowest order of 
|1 — p|, approximated to be. 


±2z 

V^TT^ 


g^ cosC^ky±'^)T 9 5sin(^fcj^±^) 




(34) 

in which the upper (lower) sign is designated for A' {B'), 
respectively. From the numerical result shown in Fig. 


the factor of p = 1.05 splits the quadratic point into two 
Dirac points located along kx axis. On further increasing 
g, a larger separation in kx occurs. On the other hand, 
for p < 1, the splitting then takes place along the ky axis. 


V. CONCLUSIONS 

The bulk entanglement spectrum of integer quantum 
Hall ground states under the checkerboard partition is 
demonstrated to be a solvable case through our formu¬ 
lation. The manifestation of topological characters, the 
filling factor v and flux <i> threading each pixel, in band 
crossings of bulk entanglement spectrum clearly shows 
the close relation between gapless entanglement spectrum 
and Chern number of underlying physical system. The 
appearance of a single Dirac point in either v = 1 and 
^ = 2tt 01 V = 2 and $ = tt as well as the quadratic point, 
or two Dirac points, in the case of = 2 and $ = 27r con¬ 
firms that the classification in physical system is carried 
over to the entanglement spectrum. Moreover, the ques¬ 
tion whether the nodal line enclosing the BZ in the case 
oi 1 / = 1 and $ = tt represents half of a Dirac point is an 
interesting one. In this paper, we only study the cases 
with N = 2 and 4 and find the spectra are all symmetric 
with respect to zero, which is consistent with the fact 
that the correlation matrix G in Eq. can be mapped 
to I — C upon complex conjugation. On the other hand, 
it is straightforward to study the cases of iV = 1,3, • • •, 
for which our preliminary results show that one band 
crosses zero but the other iV — 1 bands are gapped but 
form charge conjugate pairs. 

It deserves further study of classifying the topological 
nature of BES by using the emerging discrete symme¬ 
tries we have identified, similar to what have been done 
in classifying the gapless edge modes of topological insu¬ 
lators/superconductors. Moreover, we have only consid¬ 
ered the commensurate cases, i.e., N is positive integer, 
it is also interesting and challenging to solve the BES 
of the incommensurate cases and find out their physi¬ 
cal implication in the similar context of disorder-induced 
phase tra nsitio n for the BES of AKLT state studied in 
the works 
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